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Abstract 

For quantum mechanics of a charged particle in a classical external electromagnetic field the momen- 
tum and Hamiltonian operators are gauge dependent. For overcome this difficulty we reexamined the 
effect of a gauge transformation on Schrodinger and Dirac equations. We show that the gauge invariance 
of the operator H — ih-^ provides a way to find the energy operator from first principles. In particular, 
when the system has stationary states the energy operator can be identified without ambiguities for 
non-relativistic and relativistic quantum mechanics. Finally, we examine other approaches finding that 
in the case in which the electromagnetic field is time independent, the energy operator obtained here is 
the same as one recently proposed by Chen et al. []. 

1 Introduction 

The physical description of nature must not depend on arbitrariness arisen from mathematical formalism. 
Typically some non-physical quantities (NPQ) are introduced in order to find a solution for the equations 
which governs the dynamics of a given system. In electrodynamics a gauge transformation of the electro- 
magnetic fields is a well known example of such an arbitrariness. Classical mechanics with electromagnetic 
fields, via Lorentz force, is gauge invariant [2, 3]. However, there are apparent problems concerning the 
momentum and Hamiltonian operators of the charged particle in quantum mechanics, due to the fact that 
the expectation values of these two operators are gauge dependent [1, 4]-[21]. 

Experimental results are gauge independent, therefore it is necessary to find a gauge invariant quantity 
to represent physical energy. In the literature, there are some suggestions for this purpose. One of them 
consists in a gauge invariant energy operator which was proposed by Yang [7]. This approach frequently has 
been used in literature [8]-[16]. However, as we will show, it is possible to find counterexamples in which this 
approach does not represent the energy. A different proposal suggests to restricting the gauge function to be 
a sum of a purely spatial function plus a linear function of time [18]-[21]. We will show that this approach 
yield to contradictory results. 

On the other hand, finding a definition of gluon spin and orbital angular momentum has been a long 
standing problem [22, 23]. Recently Chen et al. proposed a gauge invariant decomposition of the total nucleon 
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angular momentum into quark and gluon constituents [24, 25]. The decomposition is based on the separation 
of the electromagnetic potential into pure gauge and gauge invariant parts, which simultaneously satisfies 
the requirement of gauge invariance and the relevant commutation relations. With this decomposition a new 
energy operator was proposed in the non-relativistic and relativistic quantum mechanics [1, 4, 5]. 

The purpose of this paper is of set forth the gauge invariance of quantum mechanics. For this we will 
show, in section 2, that the operator H — ih-j^. is gauge invariant, even for relativistic quantum mechanics. 
At the same time we will show how to derive the gauge invariant energy operator from first principles 
explicitly when the system has stationary states in section. In section 3, we will discuss other approaches. 
Additionally, we will demonstrate that for time independent electromagnetic fields, our approach is consistent 
with the decomposition of the potential into a pure gauge and a gauge invariant parts, proposed in [1, 4, 5], 
is consistent with our approach. Finally we summarize our conclusions in section 4. 

2 Gauge invariance 

Let us describe the electrodynamics using the potentials A and ^4°, the classical electrodynamics is invariant 
by the so-called gauge transformations 

A^A' = A + V X , A°^A°' = A Q -^-, (1) 

where the scalar function \ is the gauge junction. Due to the fact that classical equations of motion are given 
in terms of the fields, the transformation Eq. (1) does not change anything in the dynamics, i.e. classical 
theory is gauge invariant. 

The canonical momentum and Hamiltonian are defined by 

p = mv + qA, and H = f p - - a) 2 + qA°, (2) 

2m V c / 

respectively. Here q denotes the charge of the particle. The Hamiltonian is quantized if p is replaced 
by —ihV. After a gauge transformation Eq. (1), the expectation value of the operators in Eq. (2) are 
transformed as: 

(*'|p'l*'> = <*|p|*)+g<*|V X |*) (3) 

(VIH'IV) = (*\H\*)+q(*\^-\*) (4) 

where \^'} = e lqx ^ r ' t ^ h Notice that the expectation values of these two operators are gauge dependent. 
In order to remove the gauge dependence of the expectation value of canonical momentum, one introduces 
the gauge invariant operator, 

P = P qA. (5) 

It is straightforward to check that the expectation value of this operator is gauge invariant, since the expec- 
tation value of qA cancels the additional term in the right side of in Eq. (3). However, the commutators 
between the components of P are 

[P\P J ] = -iq{d l A 3 - d J A { ) = -iqF l K (6) 

therefore P does not satisfy the Lie algebra of canonical momentum (i.e. \P % ,Pi\ = 0), then it cannot be 
the proper momentum operator [1, 5]. 
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On the other hand, the dynamics of a quantum particle is described by its (complex) wavefunction ^(r, t), 
and its evolution is determined by the equation [ ] : 

Hy(r,t)=iH-*(r,t), (7) 

where h is the Planck constant and H is the Hamiltonian. Denoting by A — (A , A) to the external 
electromagnetic potential, for non-relativistic quantum mechanics we have the Hamiltonian 

Hs = i-(p-^A)\ q A°, (8) 



2m V c 

and so we obtain the known Schrodinger equation. Meanwhile, for relativistic quantum mechanics we have 
the Dirac Hamiltonian 

H D =a- (p-^Aj +/3m + qA°, (9) 

where a and /3 are 4x4 Dirac matrices. 

Let us consider the effect of a unitary transformation U(x) on Eq. (7), 

iMJ(x)^ = U(x)HV. (10) 

From the Leibniz's law for derivatives we have that 

d dU(y) 



Thus Eq. (10) which becomes 



U( X )HUHx)+ih?^U\ X ) 



(U(x)9). (12) 



Let us define the transformed wavefunction as ^'(r, t) = j7(x)^(r, t), hence the transformed Hamiltonian is 
given by 

H>^U( X )HUHx) + rt dU ^UHx) (13) 
in order to obtain the transformed Schrodinger (or Dirac) equation 

H'V = ih?£-, (u) 

Therefore, a unitary transformation U(x) does not affect the form of the Schrodinger (or Dirac) equation and 
the notion of probability, since ^'(r, t)\ 2 — |^(r, t)| 2 . However, the Hamiltonian has been changed by Eq. 
(13), which means that it is not measurable and cannot be seen as an energy operator [17]. The particular 
choice of U(x) = e lqx ^ r,t ^ h reproduces the gauge transformations Eq. (1). So the Hamiltonian is modified 
by the gauge transformation according to Eq. (13) and it does not constraint that their eigenvalues remain 
unchanged. Consequently, the eigenvalue problem of the Hamiltonian, in general, does not give the energy 
spectrum: The set of eigenvalues of H, denoted by a n , are different from the energy levels of the system. 
This set is defined as by the equation 

H(A°, A)u n = a n u n , (15) 
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where u n are the eigenf unctions of H(A°, A). Notice that {a n } in general are not the energy. In order to 
clarify the ideas above enunciated, consider the equation of eigenvalues of H' , 



H'(A a ', A')vj(x,t) = fyvjfat). 



(16) 



The sets {ak}kei y {0j}jei'-> are n °t necessary the same and so the index sets I and I' may be different 
(e.g., a spectrum may be discrete and the other not). Then, omitting the functional dependencies, 



and so 



H'{A ',A')vj 



H + ihuHx) 



UH(A°,A)U^ +ih^U^ 



PjVj, 



dU(x) 
dt 



Now, if we define fj(x,t) = W (x)vj(x, t) and remarking that U = e tqx ( T ' t " fi , one has 



H-q 



dt 



(17) 



(18) 



(19) 



Note that the eigenvalues (3 in Eqs. (16) and (19) are the same, thus H' and H — g ^*^*'*) are equivalent 
operators 1 . However, H and H' are not equivalent operators, thus they do not have the same eigenvalues. We 
conclude that in general the Hamiltonian does not represent the energy. On the other hand, the procedure 
performed to obtain Eq. (14) from Eq. (10) ensures that the operator H — ih^ is equivalent to H' — ih^. 
This equivalence implies that the solutions of the Schrodinger (or Dirac) equation for H and H' are the 
same. In this sense, the quantum mechanics of a charged particle in a classical external electromagnetic field 
is gauge invariant. 

So far we had demonstrated the gauge invariance of quantum mechanics. Now some questions rest: If the 
Hamiltonian does not represent the energy, which operator does? Can a time-dependent gauge transformation 
change the energy spectrum? 

First, let us review the recipe to obtain the energy spectrum of a particular system. Suppose that Eq. 
(7) allows the separation of variables as 



* fe (M) = e- Mfct/ Vfc(r), 
thus Schrodinger (or Dirac) equation is read as the eigenvalues equation of H(A°, A) 

H(A°,A)i; k (r)=a k ij k (r). 



(20) 



(21) 



These V'fe( r ) correspond to the stationary states of H. After applying a transformation Eq. (1), the trans- 
formed wavefunction is 

* fex (r,<) = e^ r <*% fc (r,t) = e^ r ^~ a ^Mr) (22) 

and the transformed Hamiltonian is H' (A' , A ') = H(A° — dtX,A + Vx)- Notice that the transformed 
wavefunction is also stationary, because | *f> kx (r, t) \ = \ip k (r)\. If we substitute these transformed wavefunction 



1 A linear operator C in a Hilbert space IK and a linear operator C in a Hilbert space W are called equivalent if there exist 
an isomorphism U of J£ onto Jt' (or an automorphism U of M if 5£ = W) such that C = U~ 1 C'U [27]. 
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and Hamiltonian in the Schrodinger (or Dirac) equation Eq. (7), it is easy to check that the time- independent 
eigenvalue equation associated is 



H A 



dx 

at'' 



dx 
dt 



We see that the left-hand side of the above equation is equal to left-hand side of the Eq. (21) 



H A 



dt'' 



• 9 § = *(^A), 



(23) 



(24) 



and therefore Eqs. (23) and (21) have the same solutions. Therefore, the operator obtained from the 
separation of variables is gauge invariant and represents the energy. We conclude that in the particular case 
in which the Schrodinger (or Dirac) equation allows a separation of variables as in Eq. (20), the energies are 
given by the eigenvalues a% of the Hamiltonian whose eigenf unctions are a product of an spatial function by 
e _M * fc *. Additionally, if another Hamiltonian (H 1 ) is related with the separable one (H) by means a gauge 
transformation, the eigenvalue equation of H' must be modified by introducing the term d/dt, in order to 
perform the separation of variables. Hence this new eigenvalue equation in terms of H' (Eq. (23)) provides 
the actual energies. 



3 Discussion about other interpretations 

As we mentioned in the Introduction, there are other attempts to resolve the gauge invariance of quantum 
mechanics, though, as we will show in this section, some of these approaches lead to wrong conclusions, 
whereas thee most recent proposal by Chen et al [1, 4, 5, 24, 25] is consistent with our present approach. 

3.1 Yang's energy-operator 

First we will consider the Yang's energy-operator, then we will discuss the restrictions on the gauge function 
[18] and we will finish by explaining the approach of Refs. [1, 4, 5, 24, 25]. 

The fact that the Hamiltonian is not a gauge invariant quantity, and so a NPQ, strikes against the 
common meaning of the Hamiltonian as the energy. A different proposal is given by Yang and it has been 
used in literature [7, 8, 9, 10, 11, 12, 13, 14, 15, 16]: in order define a gauge invariant operator representing 
the energy is introduced the Yang's energy- operator y(A), 

y(A) = H(A°,A)- qA°. (25) 

A simple application of Eq. (13) shows its invariance. However, a typical system, as the Hydrogen atom 
denies its plausibility. For instance, consider the n-th level of an Hydrogen atom 

*--l{<gv)h- (26) 

This standard result is obtained by using the Hamiltonian of the form H(A°, A — 0), with an electrostatic 
potential A (r) = — e/47reo|r|. This Hamiltonian does not depend on time, and so the solutions of its 
time-dependent Schrodinger equation may be calculated by using a separation of variables 

*(r,t)= e~ tE " /ht ^(r). (27) 
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For this case, the eigenvalue equation of the Yang's energy-operator is: 



y(A = 0)V(r, t) = [H(A°, A = 0) - eA>(r, t) = ^(v) = e^(v, t). 

Zm 



(28) 



Being e the eigenvalues of y(A). The above expression is obviously a free-particle Hamiltonian, it means 
that the spectrum is continuous!. Notice that Eq. (25) is equal to the free-particle Hamiltonian for any 
H = p 2 /2m + eA°(x). 

Now, if we start by taking the temporal gauge for the Hydrogen atom, we have that 



A a = 0, and A = 



et 



thus, the eigenvalue equation of y is 



y(A)^(r, t) = [H(A° = 0, A)ty(r, t) = (p eA) % (r) = e^(r, t). 

Zm 



(29) 



(30) 



Now if we multiply both sides of Eq. (30) by U(x) = e lex ^ r ' t '' ti , with %(r, t) — et/r and taking into account 
the identity 



c « ex (r,o/fi (P-e^) = [p-e(A + V%)] , i(;> , , 



2m 



the eigenvalues equation becomes 



P 

2m 



2 r 



2m 



Jex(r,t)/h 



(31) 



(32) 



Again we obtain a similar result as in Eq. (28). The above arguments proof the unsuitability of the energy- 
operator to describe properly the energy. Specifically, we proved that the operator y(A) in the Eq. (30) is 

2 

equivalent to operator and so they have the same eigenvalues. Consequently the operator y(A) in Eq. 
(25) is not the energy in general. 



3.2 Restrict the gauge function 

Another approach, proposed by Stewart [18], consists in restricting the gauge function to the form: 

X (r,t) = f(v) + g(t). (33) 

This constraint is equivalent to assume the separation of the operator: 

U( x ) = e iqf(r)/h e lq9(t)/h . (34) 

However, allowing this restriction leads to ambiguity when choosing the gauge potential. To see this, let 
us define it as a set of all possible electrodynamic potentials that may represent a given system. Eq. (33) 
divides the set il in equivalence classes (see fig. 1). According to Stewart any two Hamiltonians, H and H' , 
yields the same physical results only if they are related by a gauge transformation of the form (33). Under 
this assumption, if we take two different equivalence classes, we obtain two different physical results and 
as there is no way to choose beforehand which is the "real" equivalence class that represents the physical 
system. Then the restriction on the gauge function, Eq. (33), does not gives a satisfactory solution for the 
gauge invariance, because it does not decide which equivalence class describes properly the physical system. 
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Figure 1: Potentials space it. Each piece of this circle corresponds to a set of potentials that are related 
between themselves by Eq. (33). For each set we may chose a pair (A , A) to represent it, dividing il on 
equivalence classes (sets of potentials). (We here show few classes, but actually are infinite). 



3.3 Decomposition of the electromagnetic potential 

Finally, let us consider a decomposition of the electromagnetic potential into pure gauge and gauge invariant 
parts [1, 4, 5, 24, 25], 

A A pure + Ap/jj, s , A — A pure + Api iys , (35) 

with 

V x A pure = 0, V • A phys = 0, (36) 

V4° - 1 9A P ure v 2 4° - lVT\ 

v ^pure ~ c dt ' phys ~ £q ' v° '' I 

Under a gauge transformation Eq. (1) these two parts are transformed as follows, 



A' a i n v A' a A ' - A - — A ' - A (W\ 

"-pure "pure T V^i "-phys ~ "physj -"-pure ~ -"-pure i ^phys ^phys- \°°) 



In this approach the energy operator is H(A°, A) — qA pure = H(A p ' hys , A), where H — Hs for non-relativistic 
quantum mechanics, and H = Hd for relativistic quantum mechanics 2 . It is straightforward to prove that 
H — eA pure is gauge independent [7] . 

Next we will show that the energy operator H(A^ hys , A) is consistent with our analysis of first principles, 
at least for the particular case of a time independent electromagnetic field. In this case 

1 8A 

E(r) = -VA° - -V, B(r) = VxA. (39) 
c at 

2 It can be proved that the operator H(A® h , A.) in [ ] is the same energy operator H(A°, A) — &dt ^2"V ■ A in [24, 25], 
using that S/A° ure = -\d t A pure , one has A - fSi^-V ■ A = A° hya . 
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Without loss of generality we can choose A(r) and ^4°(r) to be also time independent, because the magnetic 
field is time independent. By means of introducing the decomposition Eq. (35), and using Eq. (36) one can 
derive 

E(r) = -VA°(r), B(r) = V x ^.(r), (40) 

on the other hand using Eq. (37), VA^ ure = and implies that A® re = const. This constant can always be 
absorbed by A° hys . Then, Eq. (7) can be written as 

H(A° phys (r), A(r))*(r, t) = »n^*(r, t), (41) 

that allows as to have variable separation of the kind of Eq. (20). The energy operator obtained their 
H(Ap h (r), A(r)), is the same as one recently proposed by Chen et al. [1, 4, 5]. 

4 Conclusions 

In this paper we have reexamined the gauge invariant of non-relativistic and relativistic quantum mechanics 
of a charged particle in a classical external electromagnetic field. We show how the gauge invariance of the 
operator H — iK-j^ provides a way to find the energy operator. In particular, when the system has stationary 
states the energy operator can be identified without ambiguities from first principles. We have dismissed 
some approaches suggested in the literature and we have found that the energy operator obtain in this paper 
(for time independent electromagnetic field) is the same one recently proposed by Chen et al. [1, 4, 5]. 
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